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The Paschen-Back effect in the Li 1 6708 A line and the 
presence of lithium in cool magnetic Ap stars* 
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ABSTRACT 

Context. A number of cool magnetic Ap stars show a prominent feature at A 6708 A. Its identification with Li i remains controversial 
due to a poor knowledge of the spectra of rare-earth elements that are strongly enhanced in peculiar stars and can potentially provide 
an alternative identification. 

Aims. We suggest to investigate the 6708 A line in Ap stars with strong magnetic fields. In these objects the magnetic broadening 
and splitting provides an additional powerful criterium for line identification, allowing to use the whole line profile instead of a mere 
coincidence of the observed and predicted wavelength. 

Methods. Due to a small separation of the Li i doublet components, their magnetic splitting pattern deviates from the one expected 
for the Zeeman effect even in relatively weak fields. We carry out detailed calculations of the transition between the Zeeman and 
Paschen-Back regimes in the magnetic splitting of the Li i line and compute polarized synthetic spectra for the range of field strength 
expected in Ap stars. Theoretical spectral synthesis is compared with the high-resolution observations of cool Ap stars HD 116114, 
HD 166473 and HD 154708, which have a mean field strength of 6.4, 8.6 and 24.5 kG, respectively, and show a strong 6708 A line. 
Results. High-resolution spectra for the 6708 A region are analysed for 17 magnetic Ap stars. The presence of the 6708 A line is 
confirmed in 9 stars and reported for the first time in 6 stars. The strength of the Li i doublet does not correlate with absorption features 
of any other element. The stars HD 75445 and HD 201601 provide an extreme example of the two objects which are dissimilar with 
respect to the 6708 A line but very close in the atmospheric parameters and abundances of other elements. We demonstrate that 
the observed profiles of the 6708 A line in the strong field stars HD 116114, HD 166473 and HD 154708 correspond rather well 
to the theoretical calculations assuming the Lii identification. Inclusion of the Paschen-Back effect improves the agreement with 
observations, especially for HD 154708. 

Conclusions. Results of our study confirm the Li i identification proposed for the 6708 A line in cool Ap stars. 
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1. Introduction 

Chemically peculiar A and B stars (also know as Ap/Bp stars) 
represent unique astrophysical laboratories for investigation of 
the hydrodynamic processes, rotation and magnetism in the up- 
per main sequence stars. Cool magnetic Ap stars (T^ $ 8000 K) 
are especially interesting due to their slow rotation, the presence 
of non-radial oscillations, strong global magnetic fields, and re- 
markable richness of the optical spectra. Abnormal strengths of 
spectral lines in these stars reveal large deviations of the stel- 
lar surface chemistry from the scaled solar abundance pattern. 
Modern abundance analyses (see a summary by Ryabchikova 
2005) typically find that the light elements are underabundant, 
the iron-peak elements are mildly enhanced or close to solar, 
whereas huge overabundances are reported for the heavy ele- 
ments, especially rare-earths. Adding complexity to this obser- 
vational picture, chemical elements are distributed inhomoge- 
neously in many Ap stars, showing both a non-uniform hori- 
zontal structure (chemical spots) and different distribution with 
height (stratification) in the stellar atmosphere. Combined with 
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a prominent Zeeman broadening and splitting of spectral lines, 
these unusual properties make the appearance of Ap-star spec- 
tra dramatically different compared to the normal stars of sim- 
ilar spectral classes. In particular, the line density is greatly in- 
creased, with the spectral regions accessible to typical modern 
CCD spectroscopy populated by thousands of unidentified lines 
of the rare-earth elements (REEs). This obviously creates signif- 
icant problems for identification and measurement of the lines of 
astrophysically interesting ions. 

The history of identification and analysis of lithium in Ap 
stars is characteristic of the difficulties encountered in the stud- 
ies of elements with a small number of absorption lines. The 
feature at the position of the resonance Li i doublet at A 6707.76 
and 6707.91 A (hereafter referred to as Lii 6708 A) was origi- 
nally reported for several Ap stars based on low-resolution, poor 
S IN observations (e.g., Wallerstein & Merchant 1965|. Modern 
observations with electronic detectors (Faraggiana et al. 1996; 
Polosukhina et al. 1999) suggested the presence of the resonance 
Li line in about a dozen stars, and also uncovered a conspicuous 
variation of this line with the rotational phase in a few objects. 
The lithium abundance exceeding the solar concentration of this 
element by 2-3 orders of magnitude is required to reproduce the 
observed equivalent width of the 6708 A line (Shavrina et al. 
2001a; Polosukhina & Shavrina 2007). Moreover, up to 6 dex 



lithium overabundance in the areas around magnetic poles is de- 
duced in the spatially-resolved Doppler imaging investigation 
of the Li spots in HP 83368 (Kochukhov et al. I2004bb and in 
HD 3980 (Drake et al. [2005). Such an unusual behaviour, com- 
bining an extreme lateral inhomogeneity and a large overabun- 
dance, is often found for the REEs in Ap stars, but is uncharac- 
teristic of the light elements. Consequently, the Li i identification 
of the 6708 A feature has been often questioned (e.g. Nesvacil et 
al. 120041 1. citing incompleteness of the currently available REE 
line lists. Indeed, an element identification hinging on a single 
line can be incorrect and resulting abundance can be significantly 
overestimated if an unknown REE line overlaps by chance with 
the expected position of the resonance lithium line. Interestingly, 
recent extension of the Cen line list has solved a similar "Li 
problem" for the s-process enriched low-mass post-AGB stars, 
where the absorption previously attributed to the redshifted Li i 
line has turned out to be due to Ce n 6708.099 A (Reyniers et al. 
[2Q02T >. 

Inconclusive situation with the identification of Li in cool Ap 
stars warrants a new look at this problem. The previous studies 
were largely concerned with measuring Li abundance and inter- 
preting the line strength variations, but have not fully explored 
the line profile information available in the high-resolution ob- 
servations of the slowly rotating Ap stars. Here we propose that 
an additional interesting constraint and, possibly, an ultimate 
confirmation of the Li identification, can be obtained through 
the study of the 6708 A line shape in stars with strong mag- 
netic field. As pointed out by Mathys & Lanz ( 1995), line iden- 
tification in the Ap stars with resolved Zeeman split lines can 
be strengthened by comparison of the observed and expected 
magnetic splitting patterns. Effectively, in stars with strong fields 
not only the central line position but also the whole line profile 
shape, determined by the number and relative strengths of the n 
and cr components, can serve as a consistency check in the cases 
where line identification is doubtful. The first practical applica- 
tion of this strategy was presented by Ryabchikova et al. (2006). 
These authors used observations of the resolved Zeeman-split 
lines in the strongly magnetic Ap star HD 144897 for a major 
reclassification of the Ndm spectrum, thereby substantially in- 
creasing the number of Nd m lines with precise wavelengths and 
atomic parameters. 

The aim of the present paper is to study the formation of 
the Li i resonance doublet in the slowly rotating cool Ap stars 
taking into account its special behaviour in a strong magnetic 
field. It is well known that for the kG-strength magnetic fields 
typically found in Ap stars, the splitting of the Li line exhibits 
substantial deviation from the normal linear separation of the 
Zeeman components due to the partial Paschen-Back effect (e.g. 
Maltby 119711 Mathys 119911 ). The importance of the Paschen- 
Back effect for the profiles of individual multiplets of iron-peak 
elements was demonstrated by Mathys (1990), Landolfi et al. 
(2001) and Stift & Leone (presentation at the Vienna CP#Ap 
Workshop, September 2007), however no study has previously 
incorporated the Paschen-Back splitting in the Li i line formation 
calculations for magnetic A-type stars. 

In Sect. [2] we address the problem of the departure of mag- 
netic splitting in the Li i 6708 A line from the simple Zeeman 
pattern. We present calculation of the partial Paschen-Back ef- 
fect in the Li doublet and demonstrate how transition from the 
weak-field Zeeman to the strong-field Paschen-Back regimes 
manifests itself in the theoretical profiles of the Li line com- 
puted for the conditions typical of the cool Ap-star atmospheres. 
Sect. [3] gives an overview of the available observations of the 



Li line in Ap stars with strong magnetic fields and presents a 
spectral synthesis modelling of the 6708 A feature in Ap stars 
HD 116114, HD 166473 and HD 154708, which turn out to be 
the most suitable objects for the purpose of our investigation. 
Finally, conclusions of the study are summarised in Sect. [4] 

2. The Paschen-Back effect for the 2 S- 2 P transition 

In the usual situation of the Zeeman effect, the splitting pro- 
duced by the magnetic field is much smaller compared to the 
energy separation between different /-levels of a multiplet. In 
this case the influence of the magnetic field can be considered as 
a small perturbation to the fine-structure Hamiltonian describing 
the system in the absence of the field. However, when the mag- 
netic splitting becomes comparable to the fine-structure separa- 
tion, the perturbation theory of the Zeeman effect is no longer 
valid. In this, so-called Paschen-Back regime, the Hamiltonians 
corresponding to the spin-orbit interaction and to the magnetic 
field effects should be diagonalized simultaneously. A compre- 
hensive derivation of the theory of the Paschen-Back effect in 
atomic systems is given by Landi DeglTnnocenti & Landolfi 
(2004, § 3.4) and will not be repeated here. It suffices to say 
that the expressions for the splitting and strengths of the n and 
cr components can be derived under the assumption of L-S cou- 
pling, but evaluation of these expressions has to be performed 
numerically for all but the simplest terms. 

The Li i doublet is formed by the transition between the 
2 S and 2 P levels. The upper level exhibits the fine-structure 
splitting into the levels with / = 1/2 and 3/2, separated by 
AEj = 0.337 cm -1 . The behaviour of the energy levels as a 
function of the magnetic field strength can be characterized by 
the quantity 



Anm^cB 

where B is the magnetic field strength, B is the multiplet separa- 
tion constant (B = AEj/3 = 0.112 crrT 1 for the 2 P term in Lii) 
and other constants have their usual meaning. In the regime of 
the Zeeman effect (w<s; 1), the magnetic splitting increases lin- 
early with the field strength. For u> ~ 1, which for the Li 1 6708 A 
doublet corresponds to B ~ 1 kG, the eigenvalues start to per- 
turb each other, and linearity of the splitting pattern is lost. This 
regime is often called the incomplete or partial Paschen-Back ef- 
fect. Finally, for a> » 1 (the complete Paschen-Back effect) the 
spin-orbit interaction can be treated as a small perturbation in 
comparison with the magnetic splitting. For the Li i 6708 A line 
the complete Paschen-Back triplet splitting pattern will emerge 
for B ^ 12 kG. Each component of this normal Zeeman pattern 
is split according to the multiplicity of the electronic levels. 

These considerations show that substantial deviations from 
the linear magnetic splitting characteristic of the simple Zeeman 
effect are expected to occur in the Li i 6708 A line for the mag- 
netic field strengths of a few kG. Such fields are common at 
the surfaces of magnetic Ap stars and in the magnetic regions 
on the Sun and active late-type stars. Thus, an investigation of 
the Li i resonance doublet requires detailed Paschen-Back treat- 
ment of the magnetic splitting for essentially all types of mag- 
netic stars. However, since the experimental work of Paschen 
and Back (Back [T9T2l . little attention has been paid in the con- 
text of astrophysical studies to the unusual response of the Li 
line to a strong magnetic field. Based on the quantum mechanical 
calculations, Darwin (1 19271 1 and Darwin ( 1928 ) presented rules 
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Fig. 1. The strength and separation of the magnetic compo- 
nents of the 2 S- 2 P multiplet as a function of the magnetic field 
strength. The relative energy AE (vertical axis) is given in units 
of the multiplet splitting constant B. Both axes are divided by 
to + 1 . Dashed lines show splitting for the linear Zeeman effect. 
Shaded curves correspond to the magnetic splitting calculated 
taking into account the Paschen-Back effect. The width of the 
curves is proportional to the strength of the respective n (darker 
curves) and cr (lighter curves) components. The vertical bars 
show the magnetic field strength (in kG) corresponding to the 
fine-structure splitting in the Na D lines (upper row) and the Li i 
resonance doublet (lower row). 



for forming chains of non-linear equations from which the pa- 
rameters of the line splitting in the partial Paschen-Back regime 
could be derived. In a series of papers, published between 1930 
and 1941 (Green [T9"30"lll94lt ), Green and collaborators applied 
the Darwin's formalism to the interpretation of the laboratory 
measurements of the emission spectra of various ions in strong 
magnetic fields. The third paper in the series (Green & Loring 
1936) examined the strong field behaviour of the 2 S- 2 P mul- 
tiplet of Li i. Using the results of Darwin's calculations, Traub 
d 19681 1, Engvold et al. d 19701 1 and Maltby ( 119711 ) studied forma- 
tion of the resonance Li i line in sunspots. They have carried out 
the line profile computations for the Li i doublet under the as- 
sumption of the Milne-Eddington atmosphere. Mathys (1991) 
presented another theoretical calculation of the Paschen-Back 
splitting in the lithium line. Mathys (1990) and Nesvacil et al. 
(2004 ) stressed the importance of taking this effect into account 
in the analyses of the Li feature in Ap stars. However, several 
recent studies of the Li i line in magnetic stars, even those sup- 
posedly presenting refined modelling of the Li transition in the 
magnetic fields exceeding 3-5 kG (e.g. Polosukhina & Shavrina 
2007, Leone 2007), have failed to include the Paschen-Back ef- 
fect. 

The goal of our study is to carry out accurate line profile 
calculations for the Lii doublet in strong magnetic fields. Two 
approaches were used to evaluate the separation and relative 
strengths of the n and cr components of the Li line formed in 
the Paschen-Back regime. First, we have solved the system of 
equations for the 2 S- 2 P multiplet following the formalism de- 
veloped by Darwin d 19271 ). Then, we have performed equiva- 
lent calculations with the help of a computer program provided 



by E. Landi DeglTnnocenti, who implemented more advanced 
mathematical methods of the Racah algebra described in Landi 
DeglTnnocenti & Landolfi (2004). Results of the two alternative 
computations are in excellent agreement. In Fig. Q] we illustrate 
the separation and strength of the 2 S- 2 P multiplet components 
as a function of to lip) + 1). In this scaled format the figure char- 
acterizes behaviour of all astrophysically important 2 S- 2 P tran- 
sitions (the Li i resonance doublet, Na D, Ca H and K, Mg n res- 
onance lines). The vertical bars in Fig. Q] indicate the magnetic 
field strength corresponding to a given splitting pattern for the 
Na D lines and the Li i doublet. It is clear that deviations from 
the linear Zeeman splitting for the latter feature can no longer 
be neglected for B ^ 2 kG, while significant asymmetry in the 
intensity of the red and blue n components appears already at 
B Z 1 kG. 

We have incorporated results of the Paschen-Back effect 
modelling in the polarized radiative transfer calculations of the 
Li i 6708 A line profile. The total oscillator strength of the Li 
multiplet is adopted from Yan et al. 0998). The presence of a 
non-zero nuclear spin implies that the lithium line exhibits hy- 
perfine splitting. This effect should be treated simultaneously 
with the magnetic splitting only for the magnetic fields below 
1 kG (Traub 1968 ). In this study we investigate the Li line in 
stars with much stronger fields, therefore a self-consistent treat- 
ment of the hyperfine and magnetic splitting is unnecessary. In 
fact, the separation of the hyperfine components is so small that 
it can be neglected for all calculations presented here. The wave- 
lengths and relative strengths of the 6 Li and 7 Li line components 
are taken from Smith et al. (1998 ). Their line list is consistent 
with the very accurate theoretical calculations (Sansonetti et al. 
119951 ) and experimental results (Volz & Schmoranzera l 1 996l . 

2.1. Local profiles 

To avoid complications due to integration of the spectra over un- 
resolved stellar surface with a complex magnetic geometry, we 
first study the Li i line profile for a solar-like case. The local disk- 
centre theoretical Stokes / and V spectra of the Li doublet are 
shown in Fig. [2] These calculations, based on the polarized ra- 
diative transfer code SYNTHMAG (Kochukhov 2007a), employ the 
ATLAS9 model atmosphere (Kurucz [19931 with T eS = 8000 K, 
logg = 4.0 and Li abundance log(/Y L i/Mot) = -8.00. This 3 dex 
enhancement of the Li concentration relative to the solar pho- 
tospheric Li abundance (log^u/Mot) = -10.99, Asplund et al. 
2005 ) is representative of the abundance of this element inferred 
for the Li-rich Ap stars. The 6 Li/ 7 Li isotopic ratio is fixed at the 
solar system value of 0.07 (Anders & Grevesse 1989). For the 
purpose of comparison with the standard treatment of the Li i 
line used in all previous studies of this feature in Ap stars, Fig. [2] 
also shows line profile calculations neglecting the Paschen-Back 
effect. Polarized line profile synthesis shows that in the magnetic 
fields stronger than « 3 kG an accurate modelling of the Li i dou- 
blet is impossible without considering the Paschen-Back effect. 
For the magnetic fields stronger than about 10 kG the line profile 
shape predicted by the Zeeman regime calculations is entirely 
different from the correct Paschen-Back treatment. 

2.2. Disk-integrated profiles 

A large variety of the global magnetic field geometries can be 
considered for Ap stars. An assumption of the low-order multi- 
polar topology is often used to fit the curves of the magnetic ob- 
servables (e.g. Landstreet & Mathys 2000). However, such mod- 
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Fig. 2. Theoretical disk-centre Stokes / and V profiles of the Li i 6708 A doublet computed with the SYNTHMAG code. The three 
panels show the spectra for different angles y between the magnetic field vector and the line of sight. Calculations for different 
field strength (B- 1, 3, 5, 10, 15, 20 and 30 kG) are shifted upward for display purposes. The Stokes V profiles are also shifted 
upward by 1 .07 relative to the corresponding Stokes / spectra. The two sets of calculations are shown: synthetic spectra including 
the Paschen-Back effect (thick curves) and the calculations treating the Li line splitting in the Zeeman regime (thin curves). 



els cannot reproduce high-resolution observations of Ap stars in 
all four Stokes parameters (Bagnulo et al. 2001; Kochukhov et 
al. 2004a) and typically fail to achieve a good fit to the Zeeman 
split lines (Bagnulo et al. 2003). Furthermore, for the slowly 
rotating cool Ap stars considered here no phase-resolved high- 
resolution observations in the Stokes / or other Stokes parame- 
ters are available. In this situation multipolar models are severely 
underconstrained. 

An alternative to using purely theoretical and often inad- 
equate multipolar geometries is to minimize the number of 
free parameters describing the magnetic field structure and con- 
sider only those quantities which can be directly constrained 
by the available observational data. From the splitting of the 
Zeeman resolved lines one can determine the mean magnetic 
field strength, while the relative intensities of the tt and cr com- 
ponents contain information about the field orientation. Thus, for 



the purpose of calculation of the disk-integrated Stokes / spectra 
it is convenient to represent the field structure with a homoge- 
neous distribution over the stellar surface. This simple magnetic 
topology is characterized by a single value of the field strength 
and a single inclination with respect to the stellar surface. Here 
we parameterize the field with a combination of the radial and 
azimuthal field components. Note that the latter is tangential to 
the stellar surface and is always perpendicular to the observer's 
line of sight. The line of sight symmetry of such homogeneous 
magnetic field geometry allows one to obtain disk-integrated in- 
tensity spectra by combining the local SYNTHMAG profiles, eval- 
uated at 7 positions between the disk-centre and the stellar limb 
(Kochukhov 2007a). This simplified treatment of the magnetic 
field topology of slowly rotating Ap stars is usually adequate for 
the purpose of abundance analysis and the line profile fitting (e.g. 
Kochukhov et al. l2002al Nielsen & Wahlgren|2002l Kochukhov 
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Fig. 3. Theoretical disk-integrated Stokes / profiles of the Lii 
6708 A doublet. Calculations for different field strength ((B) = 1, 
3, 5, 10, 15, 20 and 30 kG) are shifted upward for display pur- 
poses. The spectra computed for different magnetic field geome- 
tries are shown with the thin solid line (homogeneous radial 
field), dashed line (homogeneous azimuthal field) and the thick 
solid line (centred dipole field aligned with the line of sight). 



2003, Leone et al. l2003l Shavrina et al. 20061 Rya bchikov a et al. 
[20041 [20061 . In fact, a recent study (Kochukhov |2007bJ) of the 
triplet lines in the spectra of strongly magnetic Ap stars revealed 
that the narrow cr components in a few such objects indicate a 
much smaller surface scatter of the field strength compared to 
the predictions of low-order multipolar magnetic models. This 
provides an additional justification for the simplified single field 
strength value approach to synthesizing the Stokes / spectra of 
magnetic A stars. 

Eventually, the homogeneous field assumption should fail 
when a large range of the magnetic field strength variation over 
the visible stellar surface is present, inducing a broadening of 
the cr components. The only known observational example of 
this situation in a cool Ap star is found for HD 154708 - a star 
with an extremely strong magnetic field (Hubrig et al. 2005). 
Kochukhov (2007a) has shown that a major improvement in the 
fit to the triplet profile of the Yn 5087 A line in this star can 
be achieved with a dipolar field aligned with the line of sight. 
The Stokes / spectra corresponding to this magnetic geometry 



can be modelled with the help of a simple modification of the 
disk integration procedure applied to the local profiles produced 
with SYNTHMAG. Instead of a constant field strength and field in- 
clination values in 7 concentric rings, we can adopt parameters 
appropriate for the aligned dipole. 

In this section we present calculation of the Li i line inten- 
sity for 7 values of the mean field modulus, {B), ranging from 1 
to 30 kG. The model atmosphere, atomic parameters and abun- 
dances are identical to those adopted previously. Calculations 
are carried out for the extreme cases of the homogeneous field 
structure (radial field and azimuthal field) and for the dipolar ge- 
ometry aligned with the line of sight. In the latter case the polar 
field strength of the dipole is adjusted to yield a given value of 
(B). The stellar rotation and the instrumental broadening are ne- 
glected. 

The resulting profiles of the Li line are illustrated in Fig. [3] 
(only magnetic synthesis accounting for the partial Paschen- 
Back effect is shown). It is evident that the morphology of the 
lithium line profile is not very different from the local case con- 
sidered in Sect. 12.11 although variation of the field inclination 
with respect to the line of sight and the surface changes of the 
field strength for the dipolar geometry introduce some smooth- 
ing of the calculated spectra. 

One can notice that the surface variation of the magnetic field 
strength in the dipolar model also leads to additional differential 
broadening of the cr components when those become resolved. 
But for all the three field structures considered, the characteristic 
simple triplet splitting pattern emerges in the 6708 A line for the 
mean field modulus above a; 12 kG. This behaviour can be a key 
factor to support or refute the Li line identification, provided that 
an absorption in the 6708 A region can be detected for an Ap star 
with a strong enough magnetic field. 



3. The Li line in magnetic Ap stars 

To clarify the question of the presence of Li in magnetic Ap stars 
and to choose the best objects for application of the theoretical 
framework outlined in Sect. [2] we examined the high-resolution 
spectra of 17 slowly rotating late A magnetic peculiar stars. 
All observations analysed in the present study were obtained 
with the UVES spectrograph at the ESO VLT, using either the 
dichroic setup (390+580 nm) with a 0.5" slit or the 600 nm red 
arm setting with an image slicer. The resolving power reaches 
80000 for the slit spectra and 115 000 for the observations ob- 
tained with the image slicer, respectively. Detailed description 
of the reduction procedures employed for both types of the 
UVES data can be found in Kochukhov et al. d2"U56l I2007ab . 
All observations in the 600 nm red-arm mode were acquired 
in the context of the programs devoted to the time-resolved 
monitoring of pulsations in rapidly oscillating Ap (roAp) stars. 
These spectral time series, containing typically 60-150 individ- 
ual exposures, were combined to yield S /N> 300 for most of 
the targets. Comparable noise levels were also achieved for the 
UVES slit spectra. One-dimensional extracted and merged spec- 
tra were rectified in the 6680-6730 and 6130-6170 A regions 
using low-order polynomials. The first spectral interval is em- 
ployed in the analysis of the Lii doublet; the second segment 
is useful for the determination of the average magnetic field 
strength from the separation of the resolved Zeeman components 
of the magnetically-sensitive Fe n line at A 6149 A (Mathys et al. 
[19971) . 

The sample of 17 objects studied here includes 9 sharp- 
lined Ap stars in which the presence of the 6708 A feature 
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Fig. 4. Spectra of magnetic Ap stars in the 6708 A region. Stars 
are arranged from top to bottom in the order of increasing mag- 
netic field strength. The vertical lines show position of the Li i 
doublet (dashed lines), as well as positions of some prominent 
rare-earth lines, Cen 6706.051, Prin 6706.703, Smn 6707.473, 
and Ce n 6708.099 A (dotted lines). 



was previously reported in the literature (Faraggiana et al. 
[T9961 Polosukhina et al. [19991 Kochukhov 120031 . The three 
broad-lined stars with a strong variable 6708 A line (HD 3980, 
HD 60435 and HD 83368) are not considered in our study be- 
cause their relatively large projected rotational velocities and the 
prominent horizontal inhomogeneities prohibit an accurate in- 
vestigation of the subtle magnetic broadening and splitting ef- 
fects. 

A collection of the 16 spectra of cool Ap stars in the Li re- 
gion is presented in Fig. [4] The stars are arranged according to 
the mean magnetic field strength, which is taken from the liter- 
ature or determined using the Fen 6149 A line. The spectra are 
shifted to the laboratory frame using 3-5 lines of the iron-peak 
and rare-earth elements in the 6690-6720 A region. All stars 
except HD 217522 clearly show an absorption at the expected 
position of the Li i 6708 A line. The presence of this feature in 
the spectra of HD9289, HD 116114, HD 122970, HD 154708, 
HD 193756 and HD 216018 is reported here for the first time. 

From the comparison of the 6708 A line with the strength of 
the nearby lines of Nd m, Pr m, Ce n, Gd n and Sm n one can con- 
clude that the former line does not exhibit an obvious correlation 
with the absorption due to any known REE features, thus mak- 
ing an alternative REE identification of the 6708 A line unlikely. 
The cases of HD 217522 and HD 75445, in which no Li line is 
detected, provide a particularly convincing illustration of the de- 
viating behaviour of the 6708 A line. The former star is one of 
the coolest magnetic Ap stars known. It shows strong lines of 
the REEs in several ionization stages. All four rare-earth lines 
in the vicinity of the Li i feature are prominent in the spectrum 
of HD 217522, yet no absorption at the position of the lithium 
resonance doublet is observed. 

The Ap star HD 75445 appears to be very similar to the 
bright roAp star HD 201601 (y Equ). Comparison of the stel- 
lar spectra in the 6140-6162 A region (Fig. [5] upper panel) 
shows little difference apart from the marginal effect caused by 
a slightly stronger magnetic field in y Equ. The abundance anal- 
yses by Ryabchikova et al. (|2002 2004) show that the two stars 
closely resemble each other in the atmospheric parameters and 
chemical composition. Nevertheless, the two stars exhibit re- 
markable discrepancy in the Li i region: y Equ shows a strong 
line at A 6708 A, while this line is absent entirely in HD 75445 
(Fig. E] lower panel). This finding again supports the view that 
the 6708 A feature cannot be produced by an unidentified line of 
a common REE ion. 

The sequence of spectra presented in Fig. |4] highlights the 
difficulty of quantitative analysis of the Li region. The Li line 
is often blended by Smn 6707.473 and Cen 6708.099 A, most 
prominently in HD 101065. In stars with stronger fields, the 
wings of the Li i line and the resolved Zeeman components of 
Prm 6706.703 A start to overlap, further complicating the anal- 
ysis. 

It is not the intention of this paper to perform detailed mod- 
elling of the 6708 A region and to present accurate Li abun- 
dance determination for every star in Fig. |4] Instead we wish 
to answer the question if theoretically predicted Paschen-Back 
splitting signature of the Li line can be observed, thus provid- 
ing a strong argument to support the identification of the 6708 A 
feature with the resonance Li i line. Previous studies of the Li 
region in Ap stars (Shavrina et al. 1200 lbl 120061 Polosukhina 
& Shavrina 120071 1 have succeeded in obtaining reasonable fits 
to the Li blend neglecting the Paschen-Back effect for the stars 
with fields < 6 kG. This was achieved by introducing a number 
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Fig. 5. Comparison of the spectra of HD 75445 and HD 201601 
in the 6140-6162 A (upper panel) and 6700-6720 A (lower 
panel) regions. The two stars have very similar abundances and 
the line absorption, except for the 6708 A line, which is strong 
in HD 201601 but absent in HD 75445. 



of additional fitting parameters (different types of broadening, 
field strength, field inclination were adjusted on the line-by-line 
basis, additional low-quality predicted REE lines were used). 
This means that, although the Paschen-Back effect could be im- 
portant for these stars, its observational signature is ambiguous. 
Therefore, in the context of our investigation, we will study in 
detail only a few strong-field stars, showing a particularly dis- 
tinct Li line profile shape or exhibiting resolved components of 
the Li blend. 

Magnetic splitting of the resonance Li line is not observed in 
any star except, possibly, in HD 154708 ((B) = 24.5 kG). This is 
in a qualitative agreement with the theoretical predictions of the 
behaviour of the Li line in magnetic field. In Sect. [2] we found 
that even in the most favourable cases of negligible rotation, re- 
solved components of the Li i doublet should appear only for the 
fields exceeding w 12 kG. The two other stars, distinguished by 
the strong magnetic fields and by the appearance of their spec- 
tra in the Li region, are HD 1 161 14 and HD 166473. These ob- 
jects have the third and the second strongest fields, respectively, 
among the Ap stars analysed here and, at the same time, show a 
moderately strong Lii feature. Furthermore, in HD 116114 this 
line is free from interfering absorption of the REEs. Thus, we 
choose HD 1 161 14, HD 166473 and HD 154708 for the detailed 



spectrum synthesis analysis based on the Paschen-Back calcula- 
tions of the magnetic splitting in the Li i 6708 A line. 

Our magnetic spectrum synthesis is carried out for the two 
types of the magnetic field geometries discussed above (homo- 
geneous and pole-on dipolar). The relevant magnetic field pa- 
rameters and v e sin i are adjusted with the help of lines with a 
simple and distinct Zeeman patterns (triplets). We do not attempt 
to tune these parameters further in order to improve the Li line 
description. Instead, the observed spectra are fitted by varying 
only the element abundances and macroturbulent broadening. 
The latter is required to account for the broad line profiles of 
some ions in cool Ap stars (Kochukhov 2003 ; Ryabchikova et al. 
120071) . The reality of turbulent broadening is also supported by 
the analysis of line profile variations in roAp stars (Kochukhov 
et al. l2007ab . On the other hand, the lack of diagnostic lines and 
the absence of time-resolved spectra does not allow us to con- 
sider chemical stratification and spotted distribution of Li. For 
these reasons, and because we adopt a simplified description of 
the magnetic field structure, our theoretical calculations will not 
give an ideal fit to observations. However, we believe that the ba- 
sic qualitative agreement between the observed and synthesized 
spectra is sufficient to verify the behaviour of the 6708 A line in 
strong magnetic fields. 

3.1. HD 11 611 4 

The cool Ap star HD 1 161 14 shows resolved Zeeman split lines, 
indicating a mean field modulus, (B), of 5.9-6.0 kG and a long 
rotation period (Mathys et al. 119971 1. In a later investigation 
Landstreet & Mathys (2000) adopted a 27.6 d rotation period 
and fitted the phase curves of magnetic observables with a global 
dipolar-like field, forming a small angle with the rotational axis 
of the star. At the same time, the rotational modulation of all 
magnetic observables is marginal and could still be consistent 
with a much longer rotation period. 

The atmospheric properties of HD 116114 were investigated 
by Ryabchikova et al. (2004). They determined r eff = 8000 K, 
log g =4.1, (B) =6.2 kG and found a moderately large overabun- 
dance of the iron-peak elements for this star. Compared to many 
other Ap stars in this temperature range, HD 116114 does not 
exhibit enhanced absorption due to the Ndm and Prm features. 
Consequently, the contribution of the Prm 6706.703 A line in 
the Li region is unimportant, which facilitates accurate analysis 
of the Li i line. For the field strength of HD 116114 we expect 
to see no splitting of the Lii line (see Figs. [2] and [3]), but the 
partial Paschen-Back effect is already quite pronounced: the Li 
blend is deeper, narrower and more symmetric, the red wing is 
less extended, and the central structure is absent. 

We have performed a spectrum synthesis modelling of the 
Li region of HD 116114 with the help of the SYNTHMAG code, 
and using the element abundances and model atmosphere of 
Ryabchikova et al. (2004). The resonance lithium doublet is 
treated in the same way as in Sect. [2] The atomic parameters 
of other lines in the Li region are extracted from the VALD 
database (Kupka et al. 1 19991 1. which includes the DREAM 
data for rare-earth elements (Biemont et al. I1999I I. The re- 
sulting line list is used throughout the present paper. The av- 
erage field strength and the field orientation in HD 116114 
were determined by fitting the resolved Zeeman components 
of the Fei 6705.101 A line. A good description of the ob- 
served profile of this feature in HD 116114 (Fig. [6} is achieved 
for (B) = 6.4 kG, a 40° inclination of the field vector rela- 
tive to the surface normal and the projected rotational velocity 



6705 6706 6707 6708 

Wavelength (A) 

Fig. 6. Comparison of the observed (symbols) and computed 
(curves) spectra of the cool Ap star HD 1 161 14. The solid curve 
shows calculations accounting for the partial Paschen-Back ef- 
fect in the Li 1 6708 A line, while the dashed curve illustrates the 
synthetic spectrum computed for the Zeeman treatment of the 
Li i line splitting. 

v e sin i = 3.5 km s _1 . Theoretical calculations for the Li i line us- 
ing these magnetic field parameters and v e sin i are presented 
in Fig. [6] A good agreement of observations and spectrum 
synthesis is found for log(Nu/N tot ) - -8.35 and a contribu- 
tion of the light lithium isotope similar to the solar system 
value 6 Li/ 7 Li=0.07 (Anders & Grevesse[T9S5J. A macroturbu- 
lent broadening of V m ac=7 kms~' is applied to the Lii line, 
while the Fei 6705.101 A line indicates V mac 55 2 kms~'. This 
range of V mac inferred from the lines of different elements agrees 
with the results obtained by Ryabchikova et al. (120071 1 . 

A comparison of the spectral synthesis for the Zeeman and 
partial Paschen-Back treatment of the Li i line splitting shows 
that the latter approach is more successful in reproducing the 
6708 A feature in the spectrum of HD 116114. In particular, 
the Paschen-Back calculations predict a narrow line core, which 
agrees better with the observed profile shape. At the same time, 
the Paschen-Back calculations lead to a deterioration of the fit 
in the red wing. This problem may be related to an unaccounted 
blending or to a spotted distribution of Li. Despite this uncer- 
tainty, it can be concluded that observations of HD 116114 are 
consistent with the Li i identification and confirm theoretically 
expected partial Paschen-Back behaviour of the magnetic split- 
ting in the lithium resonance doublet. 

3.2. HD 166473 

The cool Ap star HD 166473 shows an unusually large range of 
the field modulus variation with a period of about 10 yr (Mathys 
et al. I1997I I. Our observations of this star were obtained close 
to the maximum of (B) and, thus, are best-suited for the inves- 
tigation of the Paschen-Back effect in the Li i line. The lithium 
feature in HD 166473 was modelled previously by Polosukhina 
& Shavrina (120071 ). They have neglected the Paschen-Back ef- 
fect and used observations corresponding to the rotational phase 
of weaker (B). 

Gelbmann et al. (120001 1 analysed in detail the atmospheric 
properties of HD 166473. They have determined the stellar pa- 
rameters, r e ff = 7700 K and log g = 4.2, and have inferred a large 
overabundance of the REEs along with a moderate enhancement 
of the iron-peak elements. The effective temperature obtained 
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Fig. 7. The same as Fig.|6]for HD 166473. 

by Gelbmann et al. (200^ was subsequently confirmed on the 
basis of detailed analysis of the hydrogen Balmer line profiles 
(Kochukhov et al.|2002b]). 

Our theoretical calculations of the Li i line (Sect. |2]i show 
that for the field strength of HD 166473 and an insignificant line 
broadening one should expect to find a large influence of the 
partial Paschen-Back effect. In particular, for the fields inclined 
significantly with respect to the line of sight, the Li line profile 
acquires a distinct triangular shape, similar to the one observed 
in the stellar spectrum. 

Here we have modelled the Lii resonance doublet in 
HD 166473 using the model atmosphere parameters and abun- 
dances from Gelbmann et al. (2000). The model atmosphere of 
HD 166473 was calculated with the LLModels code (Shulyak et 
al. 12004). taking into account individual stellar chemical com- 
position and including the approximate treatment of the Zeeman 
effect in the line opacity (Kochukhov et al. 2005). 

The comparison of observations and the spectrum synthe- 
sis of the Prm 6706.703 A + Lii blend (Fig. suggests that 
the mean magnetic field vector should be inclined by « 70° with 
respect to the stellar surface normal for the field (B) = 8.6 kG 
to match the deep central component of the Prm line. At the 
same time, a large macroturbulent broadening, V mac = 9 km s , 
is needed to fit the line width. Due to the long rotation period of 
HD 166473, the rotational Doppler broadening can be neglected. 
The fit illustrated in Fig. [7] was produced using log(Afp r /W to t) = 
-7.20, log(jV L i/Wtot) = -8.05 and adopting Gelbmann et al. 
(2000) abundances for other elements. The oscillator strength 
of the Prm line, \oggf = -1.64, is from Biemont et al. (Q999). 
The inferred Pr m abundance is consistent with Gelbmann et al. 
d2000b estimate of log(jV Prln/ W t ot) = -7.60 + 0.38. 

Theoretical calculations including the Paschen-Back effect 
lead to a somewhat better description of the observations com- 
pared to the synthesis based on the linear Zeeman splitting 
of the Li blend components. The latter calculation is unable 
to match the narrow triangular core of the Lii line even for 
highly inclined fields. However, significant missing opacity in 
the 6706.9-6707.5 A interval and the large broadening compli- 
cates the assessment of the partial Paschen-Back effect signature 
in HD 166473. 

3.3. HD 154708 

A very strong magnetic field, (B) — 24.5 +1.0 kG, was detected 
in Ap star HD 154708 by Hubrig et al. ( 120531) . Their T eS deter- 



mination led to contradictory results, with the El — G temper- 
ature calibration of Geneva photometry (Hauck & North 1993 ) 
giving r e ff = 6800 K, while Stromgren photometry (Napiwotzki 
et al. 19931) pointing to r e ff =7500 K. We find that the hydro- 
gen Balmer lines in HD 154708 are better reproduced with the 
latter T e s value and that Geneva photometric indices are known 
to be susceptible to anomalous line blanketing (Kochukhov et 
al. 120051 ). Hence, here we adopt T e g = 7500 K for the spectrum 
synthesis analysis of HD 154708. The model atmosphere for this 
very peculiar star is calculated with the LLModels code, assum- 
ing log g = 4.0. This model takes into account preliminary abun- 
dance estimates derived from fitting metal lines in the two stud- 
ied spectral regions and includes the effect of magnetic field on 
the metal line opacity. 

A comparison of the theoretical spectra for a single value 
of magnetic field strength with the observed profiles of triplet 
lines shows that a homogeneous magnetic topology, similar to 
the one used for HD 1 161 14 and HD 166473, is inadequate for 
HD 154708. While the n components in the spectrum of this star 
are rather narrow and can be fitted with v e sin/x:5 kms -1 , the 
cr components show an excessive broadening, indicating sub- 
stantial spread of the field modulus over the visible part of the 
stellar surface. To account for this property of the field structure 
in HD 154708, we have modified the disk integration procedure 
applied to the local profiles produced with SYNTHMAG in such a 
way that the resulting flux spectra would correspond to a dipo- 
lar field aligned with the line of sight (see Sect. 12.21 ). With this 
modification, a major improvement in the fit to triplet lines can 
be achieved for a dipolar field structure with B p = 30 kG (see 
Kochukhov I2007al) . For this geometry the disk-averaged field 
strength is (B) = 24 kG, which agrees with the direct measure- 
ments. 

Inspection of the Li region shows that the feature which 
could possibly be identified with the splitted Li i line is blended 
by the strong line of Prin at A 6706.703 A (Fig. [8]). In the 
6705.5-6709.0 A interval the spectrum of HD 154708 shows 5 
line components with the central wavelengths 6706.0, 6706.7, 
6707.4, 6707.8, and 6708.4 A. The first two features can be at- 
tributed to the groups of the blueshifted cr and n components 
of the Prm line. The other components contain significant con- 
tribution of the absorption due to the Lii line, which exhibits 
nearly complete Paschen-Back splitting for the field strength 
of HD 154708. Adopting the magnetic field geometry described 
above, we are able to identify the feature at A 6707.4 with the 
blend of the redshifted Pr in cr components and the blueshifted cr 
component of the Li i triplet. The remaining two features belong 
to the jt and to the redshifted cr components of the Li line. 

The spectrum synthesis results presented in Fig. [8] demon- 
strate a reasonable agreement of the observed profile of the 
Prm and Lii blend with the theoretical calculations includ- 
ing the Paschen-Back effect for Lii 6708 A. These spectra 
were obtained adopting the solar system ratio of the Li iso- 
topes and using an enhanced total abundance of this element, 
log(N u /N tot ) = -8.20. The REE abundances, log(AWMot) = 
-7.70, log(AWAW) = -8.20, log(AWMot) = -8.50 were 
used. The relevant REE line parameters (Biemont et al. 1999) 
were extracted from VALD. The optimal description of the ob- 
served spectrum is achieved for V mac = 7 kms -1 . 

As expected, computation for the Zeeman splitting pattern 
of the lithium doublet produces unrealistic spectrum showing 
no resemblance to the actual observations. However, even cal- 
culations with the detailed Paschen-Back treatment of the Lii 
line do not reproduce some details, in particular the narrow 
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Fig. 8. The same as Fig.|6]for HD 154708. 

6707.8 A component, in the HD 154708 spectrum. We attribute 
this problem to the simplistic magnetic field topology adopted 
in our polarized radiative transfer calculations. In reality the 
field structure probably deviates from a pure dipole, and spectra 
may be influenced by the horizontal and vertical chemical inho- 
mogeneities. Additional constraints on the magnetic and abun- 
dance distributions in HD 154708 require spectroscopic phase- 
resolved monitoring of this star. Nevertheless, even with the lim- 
ited data available in the present study, it is remarkable to find 
that the absorption features observed in the 6708 A region of 
HD 154708 coincide with the theoretically predicted positions of 
the magnetically-split components of the Li i doublet. The latter 
behaves unusually in the strong magnetic field of this star, and 
it is very unlikely that such a coincidence (both in the central 
wavelength and in the splitting pattern) could be due to an unrec- 
ognized blending. Thus, observations of HD 154708 lend strong 
support to the suggestion that the line at 6708 A in the spectra of 
cool Ap stars is due to lithium, which is strongly overabundant 
in the line-forming atmospheric region. 

4. Conclusions 

In this study we addressed the question of the unusual behaviour 
of the Li i resonance doublet in a strong magnetic field. Due to 
the small separation of the Li doublet components, its magnetic 
splitting occurs in the partial Paschen-Back regime for the field 
strengths typically found at the surfaces of magnetic stars. We 
model in detail the Paschen-Back splitting in the Li transition 
and present the first polarized radiative transfer spectrum synthe- 
sis calculations that take into account the partial Paschen-Back 
effect for this line. Comparing these calculations with the high- 
quality spectra of cool magnetic Ap stars, we are able to shed a 
new light on the long-standing problem of the Li i 6708 A line 
identification in Ap stars. 

The main findings of our investigation can be summarized as 
follows: 

• The presence of significant departures of the Li line splitting 
pattern from the one expected for the linear Zeeman effect 
in a moderately strong field is hereby confirmed. The local 
line profile calculations show that the Paschen-Back effect 
becomes noticeable in the Stokes V for B ~Z 1 kG and and 
in the Stokes / for B 3 kG, while a major modification of 
the Li i line profile shape is found for B > 5 kG. Since the 
magnetic fields of the Li -rich Ap stars are often stronger than 
this limit, these results question the validity of the previous 



attempts (e.g., Polosukhina & Shavrina 2007) to model the 
Li i line ignoring the Paschen-Back effect. 

• Calculations of the disk-integrated Stokes / spectra show 
that the triplet splitting pattern emerges in the Lii line for 
(Z?) <: 12 kG, offering a unique opportunity to verify the 
6708 A line identification in Ap stars with a very strong mag- 
netic field. 

• A survey of the available high-resolution UVES spectra of 
cool Ap stars reveals the presence of significant 6708 A ab- 
sorption in 15 stars, out of which 6 objects are identified as 
potentially Li-rich stars for the first time. 

• The intensity of the 6708 A line does not correlate with the 
overabundance of the REEs, thus invalidating the hypothe- 
sis that this line is produced by an unidentified feature of 
one of the common REE ions. The cases of HD 75445 and 
HD 217522, which do not show the Lii absorption but are 
very similar in their REE abundances to known Li-rich Ap 
stars, emphasize the difficulty of attributing the 6708 A line 
to anything but Li. 

• We undertook spectrum synthesis modelling of the Li i line 
in Ap stars HD 116114, HD 166473 and HD 154708. The 
latter object has the strongest magnetic field among all 
cool Ap stars and, thus, represents a unique natural lab- 
oratory for the study of the Paschen-Back effect in the 
resonance lithium line. We find that the synthetic spectral 
calculations which properly account for the Paschen-Back 
splitting improve the fit to the Lii line in the spectra of 
HD 116114 and HD 166473. For HD 154708, theoretically 
expected Paschen-Back triplet splitting of the Li i line is de- 
tected in the observed spectrum. Given the difficulty of treat- 
ing the unusually strong magnetic field of this star and the 
lack of observational constraints for the field geometry, the 
agreement between the predicted and observed positions and 
relative strengths of the resolved magnetic components of the 
blend containing the Li i line is deemed satisfactory. For all 
three stars we need roughly 2.6-2.9 dex overabundance of Li 
relative to the chemical composition of the solar photosphere 
to fit the observed line profiles. 

• In summary, we confirm the Li i identification of the 6708 A 
line in cool Ap stars. The puzzling diversity of the strength 
of the Li absorption in stars with similar atmospheric pa- 
rameters suggests that this element is sensitive to the time- 
dependent diffusion effects (Alecian 1998), or, possibly, is 
influenced by the chemical weather recently discovered in 
the non-magnetic chemically peculiar stars (Kochukhov et 
al. l2007b ). Interesting and diverse behaviour of Li in cool 
Ap stars challenges modern atomic diffusion theories (e.g., 
Alecian & Stift 2006) to find a satisfactory explanation of 
the accumulation of this element in the line-forming atmo- 
spheric layers. 
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Note added in proof 

In a publication, which appeared after submission of the 
present paper, Stift et al. (2008 ) presented a theoretical study of 
the partial Paschen-Back effect for a number of Fe and Cr lines 
observed in Ap stars. They have also calculated the Paschen- 
Back splitting for the Li i line, illustrating their results with a 
figure similar to Fig. 1 of the present paper. 



